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Motivated by recent experiments carried out by Esslinger's group at ETH [arXiv : 1 1 1 1 . 50 20 ^ , we 
developed the local effective field theories to describe the universal properties of the topological 
quantum critical points in non-Abelian honeycomb optical lattice. We find that the ETH experi- 
ments, which manipulate Dirac points by changing the lattice anisotropy, can be realized in adjusting 
synthetic non-Abelian gauge potentials. From the effective metric tensor, which is induced by non- 
Abelian gauge potentials, we have identified two different types of local effective field theories. In 
one effective field theory the Fermi surface is a Dirac point, the short range Hubbard- like interac- 
tions are irrelevant and the collective modes are all damped. In the other effective field theory the 
Fermi surface is a nodal line, where there emerge a dimension crossover which leads to the short 
range Hubbard- like interactions being marginal and the collective modes being undamped. 

PACS numbers: 03.75.Ss, 67.85.Lm, 37.10.Jk 



Introduction. — Much of the interest in ultracold 
atomic gases comes from their amazing tunability in ex- 
periment. A wide range of atomic physics and quantum 
optics technology provide unprecedented manipulation of 
a variety of intriguing quantum phenomena. Recently, 
based on the Berry phase [l| and its non-Abelian gener- 
alization [2], Spielman's group in NIST has successfully 
generated a synthetic external Abelian or non-Abelian 
gauge potential coupled to neutral atoms. Realization of 
non-Abehan gauge potential in quantum gases will open 
a new avenue in cold atom physics [sl, 0] , which also has 
deep connections with relativistic field theories simulated 
by optical lattice [sl, Q and cosmology Q . 

Control of ultracold atoms in a honeycomb optical lat- 
tice attracts lots of experimental [H, [ol and theoretical 
interests recently [10]. In this Letter, motivated by very 
recent experiments about the occurrence of topological 
transitions by manipulating Dirac points with a Fermi 
gas in a tunable honeycomb optical lattice [9[, we will 
provide a general strategy for studying quantum critical- 
ity of non-Abelian honeycomb optical lattice. In the sit- 
uation of no interactions, this problem has been partially 
studied in Ref. [11]. Starting from microscopic model, 
we identify topological quantum critical points (TQCPs) 
from corresponding energy spectrums by exactly solving 
the tight-binding model. At these TQCPs, we can de- 
velop local effective field theories (EFTs) to capture the 
interplay of non-Abelian gauge potential and atom-atom 
interaction. The essence of non-Abelian gauge potential 
is to provide the effective metric on the manifold of quan- 
tum vacuum, which setups a topological background for 
interaction effects. 

Most importantly, two different types of quantum crit- 
icality will emerge at different quasi- momentum points, 
which can be unified and classified clearly within EFT. 
With different quantum scaling behavior, we will argue 



the success of the tight-binding model and discuss how to 
observe experimentally these quantum critical behaviors. 
We also study collective modes on the honeycomb lattice 
with synthetic SU(2) gauge potential. Our work shows 
not only the manipulatable topological structure in ultra- 
cold Fermi gas systems with a non-Abelian honeycomb 
lattice experimentally, but also the EFT framework for 
quantum criticality of cold atom systems. 

Energy spectrum and topological quantum critical 
points. — To study Fermi gas in two-dimensional (2D) 
honeycomb optical lattice with non-Abelian gauge po- 
tentials, we use the tight-binding model 

Ho = -t ^ Ul,.cl{icj')co{j(j) + h.c, (1) 

where t is hopping amplitude, (i, j) means the nearest 
neighbors, and cl{icF)^c\{icr) {c.{ia),Co{icr)) create (an- 
nihilate) an fermion at site of •- and o-sublattice 
with spin a. Tunneling operator [/, which is a uni- 
tary operator, is related to non-Abelian gauge poten- 
tials A by Schwinger line integral along the hopping path 
7^exp(i|- J A-d\) [12|. In quasi- momentum space Eq. ([1]) 
becomes 



no = -tY^ ci{ka) 
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Co{ka) + /i.e.. (2) 



In order to make the problem analytically tractable, 
we choose the gauge potentials on each hopping direc- 
tion with the specific SU(2) form Ui = e*""^^, U2 = I and 
Us — e^^^y [11]. Here a and {3 are adjustable parameters, 
(<Jaj, cFy) are Pauli matrices in the spin-components, t/1,2,3 
are illustrated in Fig. [H Then atoms hopping around 
an elementary plaquette is just a unitary transforma- 
tion Uq = U1U2U3UIU2UI ^ which describes the intrin- 
sic symmetry of the spins for the fermions. By taking 
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FIG. 1. (Color online) Bipartite honeycomb lattice can be 
considered as comprises of two triangular lattices, denoted 
by dark blue dots and o, denoted by light yellow dots. 

TABLE I. Zero modes obtained from energy spectrum. We 
only list a half of them due to the time revers symmetry. In 
unit a~^, a is the lattice spacing. 

/3 location of zero modes in moment space 

(0, '-) I (-1 arccos (^i±|£2^) + ^, --^ arccos (^i±|£2^)) 
(_ 1 arccos (^IM) + |, arccos (^1±|^) + 

7^ (_i arccos (-i±|^) + f,-^ arccos (-i±|^)) 

(_i arccos ( -±|^) + |, -^arccos ( -±|^) + ^) 
(arccos (^HpIS), ^ arccos ( ^^"^^ )) 
(arccos (-^jEpEg^ jl arccos 
(^,7r) (arccos ( ^^^^f^ ) , ^ arccos ( ^^"^^ ) ) 

(arccos (-^IHpEg^ , jl arccos ^- 



trace of this unitary matrices, we can obtain the gauge- 
invariant Wilson loop W = Tr [/q = 2 — 4sin^ asin^ /3. 
The Wilson loop characterizes the non-Abelian magnetic 
flux through an elementary plaquette, which is non-local 
observable. We fix a = 7r/2 and vary /3 G [0,7r] to make 
expressions simple. The spectrum of Ho consists of four 
bands given by 

ei±(k) = ±t^h-2V~d, e2±(k) = ±t\jh^2Vd (3) 

with h = 3+2 cos /c+ cos (3^ d = (sin\/3/cy — cos (3 sin + 
sin^ P{1 — cos 3kx cos V^ky) and /c± = 3/2kx ± V3/2ky . 
The spectrum is particle-hole symmetric, where ei de- 
scribes low-energy bands and 62 describes the higher- 
energy bands. 

By solving ei+(k) = for k, we can obtain zero modes 
in analytic form which are listed in Table I. Obviously, 
varying /3 from to tt, we can find two types of merging 
of Dirac points, the corresponding critical points are 
f3 = 7r/3 and f3 = 27r/3. Explicitly, when f3 G (0,7r/3), 
we have Nd = 4 massless Dirac zero modes; when 
(3 G (7r/3, 27r/3), we have Nd = 8 massless Dirac zero 
modes; when P G (27r/3,7r), we have Nd = 4 massless 
Dirac zero modes. 

Varying configuration of gauge potentials, the sys- 
tem undergoes a topological quantum phase transition 



(TQPT). Fermi surface is a topologically stable singu- 
larity of the Green function in moment space. TQPTs 
are induced by changes in Fermi-surface topology, so they 
are beyond Landau's paradigm [l3|. The different phases 
in TQPTs are not classified by different symmetries, 
but characterized by different topologies of Fermi sur- 
face. Here these phases are defined by numbers of Dirac 
zero modes Nd- The TQPT is from Nd = 8 phase to 
Nd = 4 phase. Topological quantum phase transitions 
can be defined by the singularities of the ground-state 
energ y a s a function of the parameters in the Hamilto- 
nian [1^. To get a better understanding of these sin- 
gularities, we calculate the ground-state energy which is 
^(/^) = 4^ /Bz^'k(^i-W+^2-(k)) where ei,2-(k) is 
given in Eq. ([3]) and BZ stand for the first Brillouin zone. 
We numerical calculate S{l3) and its derivatives (Fig. [2]). 
By checking this numerical result, we find that the sec- 
ond order derivative is continuous, but the third one has 
two discontinuous points. From those results, we can con- 
clude that this phase transition is the third order TQPT, 
while in contrast, most conventional continuous quantum 
phase transitions are second order. 

In order to have a close look of TQPT, we can calculate 
Berry phase for each Dirac point or equivalent winding 
number as the topological charge [l5|. Examining these 
two critical points, we find this TQPT is due to merging 
of two Dirac points (see Fig. [3|). These two Dirac points 
have opposite Berry phase, one is +7r and the other is — tt. 
This type of TQPT can be interpreted as pair production 
or annihilation of Dirac modes. 
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FIG. 2. (a) The ground-state energy and (b) its first-order 
derivative, (c) second-order derivative, and (d) third-order 
derivative. We can see that the first-order derivative and 
second-order derivative of the ground-state energy are con- 
tinuous, while the third-order derivative begins to develop 
discontinuous. 

Local effective field theory. — In what follows, we shall 
study local effective field theory at the TQCP, i.e. /3 = 
7r/3. The 4x4 Hamiltonian Ho contains information 
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FIG. 3. (Color online) Counterplot of energy spectrum for 
non-Abelian gauge field configuration: (a) /3 = 7r/2, (b) 13 — 
27r/5, (c) P = 7r/3 and (d) (3 = 7r/4. From (a) to (d), ,two 
pairs of Dirac points become close, then merge and finally 
disappear. The system undergoes a TQPT from Nd — 8 
phase to Nd = 4 phase. 



about the higher-energy band 62, thus is not suitable to 
be described by a local effective field theory. At K point, 
which denote quasi-momentums of zero modes listed in 
Table L, we first find component for each bands by diago- 
nalizing 1-Lq{K). Then we can separate Hamiltonian into 
2x2 blocks, and adiabatically eliminate higher-energy 
band contribution to get the effective two bands Hamil- 
tonian. Now, we can expand Hamiltonian around the K 
point in terms oi k = K ^ q ioi \q\ <C \K\ ~ 1/a. The 
small parameter, which governs the expansion of the en- 
ergy dispersion, is, therefore, \q\a <C 1. The continuum 
limit is reached by setting lattice spacing to approach 
zero; lima^o'^o/<^ = '^efr(k) = T'^^kiCij. The square 
of the effective Hamiltonian then gives the space metric 
g^^ = -l,g^' = 0,g'^ = TikTjk. It is easy to known that 
equal-energy surface is quadratic forms, discriminant can 
be expressed as A = — (detT)^. Berry phase can be also 
obtained as sgn(detT)7r. So the significant character of 
our system is described by whether det T = or not. We 
find that T is nonsingular (detT 7^ 0) at i^i = (|-, 

point, singular (detT = 0) at K2 = (ff?"^^) point. 
From the effective Hamiltonian we can construct local 
effective field theory at low energies as 

<Seff = j dujd?k^!l[5^^uj - (4) 

The action can be viewed as zero density critical theory 



[igI , [13, which is our starting point to explore the in- 
terplay of non-Abelian gauge potential and atom-atom 
interaction. 

Quantum scaling and collective modes. — In the fol- 
lowing, we classify our local effective field theory by the 
determinant of T. There are two distinct cases: 

(I) Dirac point (detT 7^ 0): When the matrix T be- 
comes nonsingular, the Fermi surface is a point. Thus, 
the Fermi surface has codimension = 2, which is de- 
fined as spatial dimension subtract dimension of Fermi 
surface [18]. The EFT is a massless "relativistic" Dirac 
theory. Now, we should scale the quasi-momentum as 
d^k ~ s'^d^k when rescaling of the energy uo ^ suj^ and 
the dimension of the fields will be [^] = —2 in units of 
energy. Next we consider the influence of short-range 
interactions, which are mainly interactions in cold atom 
systems, i.e. Hubbard-like interactions. Fermionic renor- 
malization group (RG) analysis shows that these interac- 
tions are irrelevant in weak coupling. This result can be 
also inferred from codimension of Fermi surface dc > 1. 
As a consequence, weak short-range interactions should 
not affect too much the low energy properties. This 
lead to our conclusion that the single particle picture 
is reliable in weak coupling case. Things become com- 
plicated when interaction is strong coupling. Compar- 
ing with weak coupling case, in which we have a con- 
ventional semimetal with fermionic excitations at Dirac 
point, there might be an insulator with Neel order and 
some intermediate phases for strong coupling [17 1. 

We can construct noninteracting Green function from 
effective action in Eq. ([4]) 



a(°Hk,iu;n) = N„-^(k)]-i = ^ 



(5) 



where projector operators are defined as, Vs 
sTijkiCFj / Ck]^ and Ck = y^kiTikTjkkj. The bare density 
and spin susceptibility can be expressed as. 



(6) 



with = {I^cr). For short-ranged Hubbard-like inter- 
actions, the random phase approximation (RPA) correc- 
tion to the susceptibilities is given by x^^^{q^^ + = 
y[l — Urx]~^ where F is a diagonal matrix diag{l^ —1) 
[19| . The collective mode for short-range interaction can 
be determined by the pole of x^^^^ which is 



V 



(1 



to 



(7) 



where Vfx = \/^^~^^^^ ^ ^ ^fy ~ V ^21 + ^22 ^ ^ 
after a proper frame rotation. It is obvious that there is 
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no real solution of uo for Eq.([7j), so the collective modes 
are damped and do not propagate. Although hopping 
amplitude t and onsite interaction V are tunable in our 
case, we cannot obtain undamped collective mode at zero 
temperature in TQCP. Since these collective modes are 
damped in the wide region of parameter space, character 
stability Dirac excitation still dominates the low energy 
physics. In this case, our Dirac point has non-zero Berry 
phase +7r or — tt, this non-zero Berry phase topologically 
protects the Dirac point. They are generally robust under 
weak perturbation which is different from the following 
case. 

(II) Nodal line (detT = 0): When the matrix T 
becomes singular, the Fermi surface is a straight line. 
Thus, the Fermi surface has codimension dc = 1. We 
should scale the quasi- momentum as d^/c ~ dk\\sdk± 
when rescaling of the energy uj ^ suj^ and the dimen- 
sion of the fields would be [^] = —3/2 in units of energy. 
Next we consider the influence of short-range interac- 
tions, which are mainly interactions in cold atom sys- 
tems, i.e. Hubbard-like interactions. At the tree level, it 
is obvious that these interactions are marginal from RG 
theory of interacting fermions. So we need go further to 
the loop level to determine whether it is marginal relevant 
or marginal irrelevant. From the exact result of Yang and 
Yang [20[ , the RG flow of coupling constant stays zero to 
all orders, if the magnitude of original coupling strength 
is less than order unity. In this case, the corresponding 
energy spectrum becomes to e{k) = ±Vfkx, it is cru- 
cial to notice that our effective field theory depends on 
one direction only. We can make a dimension reduction 
and introduce the Hubbard-like interaction, then our 2D 
system can be viewed as a set of weakly coupled one di- 
mensional (ID) spinless Luttinger Liquid in ^/-direction. 

From Luttinger liquid theory [2l|, the dispersion of col- 
lective modes reads 

u;{q) = qVfVl-VyV,^ (8) 

where Vf = \/^n+^^ ^ ^ after a proper frame rota- 
tion and critical interaction strength Vc = 27rvf. Hopping 
amplitude t and onsite interaction V are all tunable pa- 
rameters in the cold atom systems. If we tune them in 
the region V > Vc, our collective mode is damped. But 
Vc is so large that not reachable in current experiments, 
thus we cannot observe this damped mode. Moreover, 
reaching such a huge Vc large then the Fermi scale, the 
system becomes solid [22]. It is worthwhile emphasizing 
here that the solution is a non-perturbative result. 

Nodal line comes from merging two Dirac points with 
opposite Berry phase, thus its effective Berry phase is 
zero. Nodal line is no longer topologically protected. Un- 
der weak perturbation, the nodal line may disappear. In 
some circumstances, nodal line disappears leaving behind 
the pairs of Fermi points. But in a more general case 
there is an alternative destiny for a nodal line: zeros can 



disappear completely so that the system becomes local 
gapped [l3| . 

Discussion and summary. — Very recently, Esslinger's 
group in ETH has successfully identified and manipulated 
Dirac points in the band structure by observing a min- 
imum band gap inside the Brillouin zone via interband 
transitions with an ultracold Fermi gas of ^^K atoms in a 
two-dimensional tunable optical lattice [sl . The amazing 
tunability of their optical lattice structure allows for in- 
dependent adjustment of the tunnehng parameters along 
the different directions simply by controlling the intensity 
of the laser beams, which seems suitable to experimen- 
tally explore the new features of our case (I) EFT. 

However, contrasted with case (I) EFT, the topological 
defect becomes unstable and the corresponding collective 
modes are undamped in our case (II) EFT. In high- Tc su- 
perconductor experiments, angle-resolved photoemission 
spectroscopy (ARPES) is a powerful tool to visualize the 
nodal line (point) defects. Although the technique of 
ARPES has not yet been realized in cold atom exper- 
iment, there is a very similar technique of momentum- 
resolved radio- frequency spectroscopy [23| , which will be 
a promising tool to detect nodal line defects in ultra- 
cold Fermi gas. To further test the interaction effect on 
collective modes, one needs to have evidence for their 
linear dispersion relation dependence on the interaction 
strength, which is well tunable with the technique of Fes- 
hbach resonance. Ultimately, to make contact with ex- 
periments involving nodal line of ultracold fermions in 
a non-Abelian optical lattice, we must take into account 
the atom-atom interactions and attempt to determine the 
critical value of interaction strength Vc . 

In summary, we construct a universal EFT to de- 
scribe Fermi gas in a non-Abelian honeycomb lattice. 
This approach enable us to unify different universality 
class emerged from TQCPs and study their collective 
modes under the interplay of non-Abelian gauge poten- 
tial and atom-atom interaction. We show that ETH 
experiments, which manipulate Dirac points by chang- 
ing the lattice anisotropy, can be realized in adjusting 
non-Abelian gauge potentials. Our results are of partic- 
ular significance both for manipulating topological struc- 
ture in ultracold Fermi gas systems with a non-Abelian 
honeycomb lattice experimentally and for further study- 
ing quantum criticality of cold atom systems within the 
framework of EFT theoretically. 
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